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In this paper we analyze the canonical forms into which 
any pure three-qubit state can be cast. The minimal forms, 
i.e. the ones with the minimal number of product states built 
from local bases, are also presented and lead to a complete 
classification of pure three-qubit states. This classification is 
related to the values of the polynomial invariants under local 
unitary transformations by a one-to-one correspondence. 

PACS Nos. 03.67.-a, 03.65.Bz 



I. INTRODUCTION 

Non-local quantum correlations or entanglement be- 
tween space-separated parties is one of the most fertile 
and thought-generating properties of quantum mechan- 
ics. Recently it has become a very useful resource for 
many of the applications in quantum information theory 
and this has led to a lot of work devoted to understanding 
how it can be quantified and manipulated. 

Bipartite pure state entanglement is almost completely 
understood, while many questions are still open for the 
mixed state case. For pure states, the Schmidt decom- 
position |l| has proven to be a very useful tool, since 
it allows to write any pure state shared by two parties 
A and B in a canonical form, where all the information 
about the non-local properties of the state is contained 
in the positive Schmidt coefficients. The non-local prop- 
erties of quantum states can be also specified by means 
of other quantities invariant under the action of local 
unitary transformations. An interesting type of these in- 
variants are given by polynomial combination of the co- 
ordinates of the state in a product basis, and the relation 
between these invariants and the Schmidt coefficients is 
well known. 

Some novel aspects, compared to the bipartite case, 
appear for entangled systems of more than two par- 
ties. In this work we study the canonical forms of three- 
qubit pure states, extending the results of bipartite sys- 
tems. First we analyze the forms proposed for gener- 
alizing the Schmidt decomposition for three-qubit pure 
states. Then, we relate one of these decompositions to 
the polynomial invariants studied in [g-^. We give a 
one-to-one correspondence between a canonical form for 
a three-qubit pure state and a complete set of polyno- 



mial invariants describing its entanglement properties. 
We also classify the different types of canonical forms 
by means of the minimal number of local bases product 
states (LBPS), i.e. the minimal number of non-local pa- 
rameters, needed for the specification of a state. For any 
three-qubit pure state we give its decomposition with the 
minimal number of LBPS and the procedure that has to 
be applied in order to build it. Finally we indicate how 
to generalize the results to systems of A^-qubits, where 
many difficulties arise. 



II. GENERALIZATION OF THE SCHMIDT 
DECOMPOSITION 

The Schmidt decomposition has been a very useful tool 
for the study of entanglement properties of bipartite sys- 
tems. For a generic bipartite pure state 1$) G C^^ ® C^^ 
it reads 



m = j2a. 



a, > 0, 



(1) 



where I = min((ii, ^2), \ii) = \i)A '9 |'*)b, being |i) or- 
thonormal vectors in each subsystem, and ai are the 
Schmidt coefficients. It would be very interesting to 
find for three-qubit pure states a canonical decomposition 
generalizing the features of the Schmidt decomposition. 
However, the trivial generalization is not possible llG] 
and it is not evident how to extend the Schmidt decom- 
position to the case of A^-party systems (A^ > 2). Indeed 
several forms have been proposed (see for instance Uu). 
In recent work we gave a generalization of the 
Schmidt decomposition for three-qubit pure states, in 
the sense that the coefficients of this decomposition carry 
all the information about the non-local properties of the 
state, and do so minimally and unambiguously, i.e. the 
decomposition is not superflous. Starting from a generic 
state shared by three parties. A, B and C, 



l*> = I]i»jfeKjfc>, 



(2) 



i,j,k 



where \iik) = \i) a®\J) B®\k) c , we look for the local bases 
that allow to write (0) with the minimal number of LBPS. 
A simple counting of parameters shows that at least five 
product states built from local bases are needed in order 
to specify a generic state belonging to C^(8)C^(8)C^. There 



are three inequivalent classes of five LBPS: the first one 
is the symmetric set 



{|000),|001>,|010),|100>,|111>}, 



(3) 



the second is weakly asymmetric and corresponds to the 
three sets of states, 



{|000),|001),|100),|110),|111)} 
{|000), 1001), 1011), |100),|111)} 
{|000),|010),|100),|101),|111)}, 



(4) 



and the third one is strongly asymmetric, and corre- 
sponds to the three sets 



{|000),|100),|101),|110),|111)} 
{|000), 1010), 1011), |110),|111)} 
{|000),|001),|110),|101),|111)}, 



(5) 



where the three sets of states of (^ are related by permu- 
tation of the parties, and the same happens for the sets 
The non-equivalence between the sets (||), (Q) and 
follows from the different degrees of orthogonality be- 
tween the five states within each set (see figure 1). In g] 
it was proved that any three-qubit state can be written 
in terms of the product states of any of the asymmetric 
sets. Let us sketch the procedure. 

Starting from a generic state (g), one introduces the 
matrices Tq and Ti with elements 



m) 



jk 



ijk- 



(6) 



A change of basis on the first qubit transforms these ma- 
trices in the following way. 






(7) 



where uf^ are the elements of a unitary matrix, while the 
effect of a change of basis in B (C) implies that each Ti is 
left (right) multiplied by a unitary matrix U^{U'~^). The 
unitary transformation on party A is chosen such that 



det(r(5) = 0. 



(8) 



There are always two solutions for this equation since (|^) 
is equivalent to 



det(ro-|-a;Ti) =0, 



(9) 



where x = -% is an unbounded complex number. Now 

"oo 

we apply two unitary matrices on parties B and C in 
order to diagonalize Tq. These operations lead to the 
matrices 



BrpljrC 



Mo = WT^U 



Ao 




A3 M 



(10) 



where Xi are real and positive, since all the phases have 
been absorbed by phase redefinitions of |0)a, |1)a, |1)b 
and |l)c ■ By means of these unitary transformations we 
have been able to write the initial state (g) in terms of 
the products states appearing in the first set of (0), i.e. 

I*) = AolOOO) + Aie*'^|100) + AajlOl) + A3IIIO) + A4IIII) 

(11) 

Equation (g) has generically two different solutions, x 
and X, so two different decompositions {Kw are possible 
for the same state |^) . By limiting the range of the phase 
factor to < 93 < TT a unique solution is found when 
< (^ < TT (see [^ for more details), and then we have a 
unique canonical form in which to cast almost any three- 
qubit pure state. For the remaining ones, when 1^ = 0, tt, 
two canonical forms exist in general; we will break this 
remaining degeneracy taking, for instance, the form with 
the smallest Ai, or, if Ai is unique, taking the form with 
the smallest Aq. It is important also to note that we have 
singled out party A in obtaining { \\/\\ ) , but we could have 
chosen any of the three parties. 

From ( |l l| ) and by applying a unitary transformation 
on the third qubit. 



|0') 



^AfTA 



(Aie''^|0)-fA2|l)) 



(12) 



it follows that any state can be written, after removing 
the phases of four of the coefficients, as, 

I*) = ?7oe*'^|000) + ?7i|001) + ?72|100) + 773IIIO) + r?4|lll), 

(13) 

with rji real and positive, which corresponds to the first 
set in (|). 

Recently, it has been shown |lj] that the symmetric 
decomposition using the set of states (||) is also possible. 
The proof is based on the fact that if a given state |^) is 
written in a basis such that the state |111) is the one that 
maximizes the overlap of |^) with any product state, i.e. 



|fiiiP=max|(*|a/37)|^ 



(14) 



the coefficients ino, ^loi f^nd foil rnust be zero (otherwise 
one could find a product state with a larger overlap). 
Therefore any state can be written as 

I*) = Kne'^lOOO) + KilOOl) + K2IOIO) + K3IIOO) + K4IIII), 

(15) 

with Ki real and positive and < < tt. Neverthe- 
less, the conditions under which the decomposition (Hq) 
is unique are not known. 

A different decomposition, which can also be thought 
as an alternative generalization of the Schmidt decompo- 
sition for three-qubit states, could be writing the state 



as a superposition of two product states, not necessarily 
orthogonal, 



|*>-a|000)+/3e'^|(^i(^2<^3) 



(16) 



with a and (3 positive real numbers. This decomposition 
is only possible when J4 j^ (see below for the definition 
of J4) Pjl^], which corresponds to the GHZ-class in p3[ , 
and it has been proved very useful for the obtention of 
the optimal GHZ distillation protocol 04]. 



III. THE SET OF POLYNOMIAL INVARIANTS 

The space of states of three qubits is C"^ ®C'^ (E)C'^ , which 
depends on sixteen real parameters (including the norm 
and the global phase). Two states, j^I'i) and 1^2), are 
equivalent, as far as their entanglement properties are 
concerned, when they can be transformed one into the 
other by local unitary transformations. Therefore the ac- 
tion of the elements of the group U{l)x SU{2) x SU{2) x 
SU{2) define orbits in the space of states, each orbit be- 
ing the equivalence class of all the states having the same 
non-local properties. Thus, and as it is well-known, the 
dimension of a generic orbit for the case of three-qubit 
pure states is ten 0], so six entanglement parameters 
should be enough to discriminate between two different 
orbits. Since the decomposition ( pT| ) is unique, it gives 
six quantities invariant under local unitaries, the five co- 
efficients Xi and the phase (p, which allow us to check 
whether two generic states belong to the same orbit, i.e. 
whether they can be connected applying local unitary 
transformations. These parameters can be thought of as 
the entanglement coordinates. An alternative, though 
two-fold degenerate, set of entanglement parameters is 
given by polynomial combinations of the coefficients tijk 
which are invariant under the group of local unitaries 
[g-g[. In this section decomposition ( [ll|) will be related 
to these polynomial invariants. 

For bipartite pure states, |$) G C^^ C'*^ a com- 
plete set of polynomial invariants, which allows to know 
whether two bipartite states have the same entanglement 
properties, is given by 



^^(Pa) = tr(/OB) I = l,...,min(di,d2), 



(17) 



where pA = trB|<i>)(<i>| and pB = tr^|<i>)((f>| are the local 
density matrices. Since the eigenvalues of these matrices 
correspond to the square of the Schmidt coefficients (|^), 
we know the relation between the polynomial invariants 
and the Schmidt decomposition M. 

As it has been mentioned above, the space of entangle- 
ment parameters of pure three-qubit states has dimension 
equal to six, so at least six linearly independent polyno- 
mial combinations of tijk invariant under local unitary 
transformations are needed in order to specify the non- 
local properties of a state, or the orbit which it belongs 



to. In 10] the six linearly independent polynomial in- 
variants of minimal degree were found. The norm is a 
trivial one, so we will not consider it and in the rest of 
the paper we will restrict ourselves to the space of nor- 
malized states, the number of non-local parameters being 
reduced to five. This implies that we have '^- Xf = 1 in 
(O). Apart from the norm, the polynomial invariants 
given in R] are 

i < /i = tr(p^) < 1 

1<I2 = tr(p|) < 1 

l<h = t^ipl) < 1 

-<l4 = tr(p^ (g) pbPab) < 1 

0</5 = ]Hdet(t,jfe)]2 <— , 



(18) 



where 



PA = trBcl*)(*l 
PB = tr AC ]*)(*] 
PC = tr AS ]*)(*] 
PAB = trc]^)(*]. 



(19) 



and Ildet(tijfc) is the hyperdeterminant of the coefficients 
tijk (iq] and corresponds to the three-tangle of Jp^. An 
equivalent set of invariants can be constructed 

1 



Ji = -{l + h-l2-h~2y/h) 



J2=-{l~h+h-h-Wh) 



J3 = -{l~h-h + h-2^h) 



J4 = 



J5 = -(3 - 3/1 - 3/2 - /3 + 4/4 - 2V/5), (20) 

which, in terms of the parameters of the decomposition 
(|ll|), are equal to 

0< Ji-]AiA4e''^-A2A3]'<i 

1 
< J2 = MO/^2 < T 

1 
< J3 = MoA^a < T 

1 
< J4 = M0M4 < T 

1 2 

-j^ <J5= Pa{Ji + P2P-. - PiPi) < ^, (21) 

where fj,i = A|. It can be proved that J4 and J5 are 
invariant under permutation of the parties, because so 



are 2/4 — Ii — I2 and I5, and Ji, J2, and J3 single out 
parties A, B and C respectively, and transform among 
themselves under party permutation. 

From the above expressions one can prove the tighter 
bounds 

0<J2 + J3 + Ji<2 
< Ji + ./3 + Ji < T 
< Ji + J2 + ^4 < T 
0< J4 + J5 < -. 



(22) 



Also the following holds 



Jl 


= 0^ 


J5 


= 


J2 


= 0^ 


•h 


= 


■h 


==0=^ 


•h 


-0 



J4 = 0^ V^i^2J3 = y. 



(23) 



From (El]) and using the normalization condition 
^j/Xi = 1, it is possible to obtain the value of the set 
of coefficients {/^i}, 



^4 



J4 + J5± 



2( Jl + Ji) 



± Ji 



: 2,3,4 

J2 + J3 



J4 



Mo" 



(24) 



where 

Aj = {Ji + J5)2 - 4(Ji + J4)(./2 + J4)(J3 + J4) > 0, 

(25) 



which implies 



J4 + J5 = ^ J4 = J5 = 0. 



(26) 



Note that the value of cos(/7 can be also found from (21) 



cos if 



fifnt + M^M^ - Jl 
2Aj A2 A3 A4 



(27) 



and thus almost all the information about the decom- 
position can be extracted from the values of the {Ji}. 
There remains however some ambiguity in these expres- 
sions, since there are two solutions for the coefficients, 
corresponding to Hq and fj,Q , and for each of them, two 
different angles, < (p < tt and (f^ — 2tt ~ (p^, com- 
ing from ( pTj ) . Part of this uncertainty is due to the two 
solutions of (||) and in fact the coefficients {fj,f,ip'^} and 
{li~ ,<p~} describe the same orbit, and the same happens 



for {fi~ ,(p~} and {fif ,(f^}. As it has been said, the so- 
lutions associated to (f are not considered because of the 
range of the angle. However the set of invariants {J,} 
(or {h}) does not determine a unique orbit, or equiva- 
lently a canonical point representing it. Two candidates 
are possible, {/i^ ,(/?^}, so there is still some ambiguity 
left. 

The five polynomial invariants (Iq) are real, and this 
means that they can not distinguish among the orbits 
associated to a given pure three-qubit state |^), with 
coefficients tijk, and to |^)*, given by t*-j,. Indeed, 



hm*) = i,mr = i,m), 



(28) 



where the second equality comes from the fact that the 
invariants are real. It is not possible, due to this ambigu- 
ity, to individuate a unique canonical state representing 
an orbit from the invariants (p^, or (20). A twelfth 
degree complex polynomial invariant, Jg, introduced by 
Grassl [|l^, solves (albeit redundantly) this problem, just 
by inspection of the sign of its imaginary part (in other 
words, the second equality of (pq) is not valid for this 
invariant). The explicit form of Grassl's invariant, using 
decomposition ( pT| ) is 

h = MoA*4(A4(l - 2(^0 + Ml)) + 2AiA2A3e-*'^)2. (29) 

The set given by (^SJ) and Iq is complete, it allows to 
check when two states belong to different orbits, and 
from their values one can obtain a unique canonical point 
representing the orbit applying ( p4[]27| ) and, in the end, 
using Iq to discriminate between the two candidates. 

This situation is quite different from what happens for 
pure states of bipartite systems. In this case, a generic 
state 1$) £ C'l (XiC^^, with coefficients i^, can be always 
tranformed into |$)* by local unitary transformations, 
as this is clear from the fact that all the Schmidt coeffi- 
cients are real. In general this is not true for three-qubit 
systems, although in some cases the state jv?) and its 
complex conjugate j^P)* are in the same orbit. This cor- 
responds to the situations when either 



I cos if ' 



I C0S(y5 



or 



COS (p'^ = COS ip 

^4 = Mr- 



(30) 



(31) 



Equivalent conditions in terms of the invariants {Ji} can 
be obtained, giving 



V J1J2J3 



\J^\ 



for the ffi'st case and 



A , = 0, 



(32) 



(33) 



for the second. Furthermore in both situations a product 
basis can be found for which all the coefficients tijk are 
real. For the states satisfying the ffist condition, this 
basis is the one that gives decomposition (|0), since we 



have 



±1, while in the second case the proof is a 



bit more tedious and it is given in the appendix A. From 
these results, then, it follows that 



ksl 



I*) - I*)* ^ W1J2J3 = -^ or Aj = ^ I*) real, 

(34) 

where a pure state belonging to C'^ (E) C'^ 1^ C^ is said to 
be real when there exists a product basis where all the 
coefficients are real. 

To summarize, five independent quantities invariant 
under local unitaries are needed to specify the non-local 
properties of a generic three-qubit pure state. The coef- 
ficients appearing in the decomposition (O) form a com- 
plete faithful and minimal set of such invariants, when 
constrained as explained after (pT|). The polynomial in- 
variants given in dlH) must be completed with /g in order 
to solve the ambiguity between the orbits associated to 
1^) and 1^)*, and from the values of these polynomial 
invariants one can build a unique canonical point rep- 
resenting the orbit. Also when |^) and |^)* are in the 
same orbit there exists a product basis where all the co- 
ordinates of 1^) are real, as it happens for pure states of 
bipartite systems. 

Let us mention finally that any real state can be writ- 
ten with real coefficients in terms of a set of six LBPS, 
adding the state |011) to (|) or to the ffi-st of (|). This is 
done by diagonalizing Tq with two orthogonal matrices. 



A. !/ = 4 

There are several subfamilies of states that allow for a 
decomposition in terms of four LBPS. 

Type 4a: This subfamily is given by the states with 
IJ.4 — in (pT|). It is easy to prove that this condition is 
equivalent to J4 = (we will take the rest of invariants 
different from zero, unless otherwise specified). Condi- 
tion (p3) is also satisfied with J5 > 0, since all the phases 
can be absorbed. 

Type 4b: States with ^2 = (^3 = 0) in (111]). The 



equivalent conditions in term of the invariants are J2 — 
J5 = (J3 = J5 = 0). Let us mention that there is an 
apparently lack of symmetry in this subfamily, but this 
is due to the fact that party A has been singled out in 
the determinations of the decomposition (|ll|) . In fact the 
analogous states with Ji = J5 = are written with four 
terms if either party B or C is singled out in (|7|-[10|). 

Type 4c: States with /ii = in (|ll|). It can be proved 
that the corresponding conditions in terms of the invari- 
ants are J1J4 -I- J1J2 + J1J3 -I- J2J3 = VJ1J2J3 = ^■ 
Again the lack of symmetry is due to the fact that party 
A is privileged in the calculation of the decomposition 
(|ll|). Analogous condition can be found interchanging 
the role of the indices 1, 2 and 3, which means that the 
minimal decompositions is obtained if one of the other 
two parties is singled out in (plO). 

Type 4d: States with kq = in (p^). It is proved in 
appendix B that the corresponding condition, apart from 
(p^), which is always satisfied when i/ < 5, is Aj — 0. 



B. !/ = 3 



IV. MINIMAL DECOMPOSITION 

We have seen that a generic three-qubit pure state can 
always be written in terms of five product states from 
any of the sets of states in (||) , (Q) or (^) . However it is 
not clear which set should be used to find the minimal 
decomposition, that is, the one with the least number 
of non-local parameters. The minimal number of LBPS 
needed to specify a state |^) will be denoted by i^(^). We 
know that in general i> — 5 but now we want to analyze 
the cases in which i/ < 5. In this section we present 
a complete classification of the three-qubit pure states 
according to this minimal number of product states. We 
also give necessary and sufficient conditions written in 
terms of the invariants {Ji} to be satisfied by the states 
of each class. The number of non-local parameters in 
each family is I'—l, since all the coefficients arc real. All 
the families satisfy condition (|3^). 



Now we move to the study of those states that can be 
expressed as a sum of three LBPS. 

Type 3a: This subfamily is given by taking ni — fj,^ — 
in ( |ll| ) . The equivalent conditions for the invariants are 
J4 = and J1J2 + J1J3 + J2J3 = VJ1J2J3 = ^• 

Type 3b: These states correspond to the case fij — 
Mfc = in (|l|), for j,k € {1,2,3} and j ^ k. These 
conditions expressed in terms of the invariants are Jj = 
Jk^Js^ 0. 



C. 1^ = 2 

The states with two product states built from local 
bases are just in two classes. 

Type 2a: J, = except Ji(J2,^3), and these are 
the states where party A(B,C) is not entangled with the 
other two parties, so there is not truly three-party entan- 
glement. 

Type 2b: Ji = except J4, they include the standard 
GHZ state. 



D. !/= 1 



Typel: Ji = 0, and these are the product states where 
there is no correlation between the parties. 



E. Summary 

All the states belonging to C^ (g) C^ (g) C^ have been 
classified in terms of the minimal number, v, of LBPS 
required to express the state, and the resulting families 
of states are shown in table |. Generically five terms are 
needed, although there are cases where v <'a. Necessary 
and sufficient conditions in terms of the set of invariants 
{Ji} are given, which can be used to recognise the sub- 
family a three-qubit pure state belongs to. Once this has 
been done, we have provided the procedure that has to 
be applied in order to find this minimal decomposition 
with product states. 



V. GENERALIZATION TO MORE PARTIES 

The decomposition (O), which generalizes the bipar- 
tite Schmidt decomposition, has been proved to be very 
fruitful for the case of three-qubit pure states, so it will 
be convenient to know the way it can be generalized to 
more parties. In this section first we will consider with 
some details the case of four-qubit systems and this will 
give us insight into the difficulties found when we try to 
extend our results. 

The procedure to be applied for the generalization of 
decomposition ( |ll| ) for pure states belonging to C'^^C'^® 
C^^C^, i.e. states |*) = ^- ■ j, ; tijki\ijkl) shared by four 
parties A, B, C and D, will be now described. First we 
define the two hypermatrices |15| 

{Ti)jki = tijki, (35) 

which means that the initial state is interpreted as 



|vI')^|O>|0o> + |l)|0i 



(36) 



where \(f>i) are, up to normalization, three-qubit pure 
states, their coordinates being given by the elements 
of the corresponding hypermatrix T^. The effect of the 
change of local bases is very similar to the one described 
for three-qubit systems: a unitary transformation on sys- 
tem A mixes the coordinates of the two \4>i), while uni- 
tary transformations on the rest of subsystems can be 
used to make zero some of their coefficients. Now we 
apply the change of local bases on A that gives 



Hdet(To) = 0, 



(37) 



found for three-qubit pure states. Since ( pT\ ) is verified, 
it is known that |0o) belongs to, at least, type 4a states, 
so we will manage to write the initial state Ivt) in terms 
of the twelve product states: 



|0000), 10100), 10101), 10110), 
|1000),|1001),|1010),|1011), 
|1100),|1101),|1110),|1111). 



(38) 



and afterwards unitary transformation on B, C and D are 
used to write the new \(J)'q) in the canonical decomposition 



A simple counting of parameters gives that the min- 
imal number of LBPS needed to specify a state |\1/) G 
C^(8)C^(8)C^(X)C^ is exactly twelve. The decomposition we 
have found depends on twenty-four non-local parameters 
but it is known that by phase redefinitions, i.e. acting 
locally with C/(l), five phases can be absorbed (generi- 
cally, for N parties A^ -I- 1 coefficients can be made real), 
so the number of non-local parameters is nineteen (in- 
cluding the norm) , as it was expected O] . 

However some problems arise in this case. Many de- 
compositions in terms of the set of states (^^ are pos- 
sible for the same state. In fact (37) is a fourth degree 
equation, so four solutions will be found and from these 
solutions four different decompositions will be derived. 
For the case of three-qubit pure state there were two 
solutions for (g|), but we managed to obtain a unique 
decomposition by limiting the range of i^. A similar rea- 
soning seems not to be trivial for this case. Furthermore 
for pure four-qubit states more inequivalent set of twelve 
product states appear, and this will difficult the analysis 
of the minimal decomposition. The generalization of de- 
composition (|T|) to A^-qubit pure states (A^ > 3) is then 
quite cumbersome. 

Finally, it has to be noted that the algorithm proposed 
in |jl2] for the decomposition dl3) can be also extended 
to higher dimensional systems. Let us mention however 
that, in any case, as the dimension of the space increases, 
the number of coefficients that can be made equal to zero 
in any of the decompositions becomes irrelevant. 



VI. CONCLUSIONS 

In this work we have studied the canonical forms of 
pure three-qubit states, extending the known results of 
bipartite systems. 

First we show the possible generalizations of the 
Schnfidt decomposition and we relate one of these de- 
compositions (O) to the polynomial invariants of P-p|. 
The six linearly independent polynomial invariants of W] 
are not able to discriminate betwee the entanglement or- 
bits associated to a state and its complex conjugate in a 
product basis. An additional polynomial invariant intro- 
duced in pTJ has to be used, and we have seen how to 
connect this complete set of polynomial invariants with 
our generalization of the Schmidt decomposition. Indeed 
it is shown how to find a canonical point in a generic 



orbit described by this complete set of invariants. Let 
us mention here that a three-qubit pure state \'^) and 
its complex conjugate l^")* give the same optimal prob- 
ability of distilling a maximally entangled state of three 
qubits, in the single-copy case [ p^ . 

We have also looked for the decomposition of any state, 
I^E"), with the minimal number, J^(^), of product states 
built from local bases. Generically this number is equal to 
five, although many exceptional states have been found 
with ly < 5. We have been able to give a complete clas- 
sification of these states by means of a set of necessary 
and sufficient conditions written in terms of the set of 
invariants (P0|). The procedure to be applied in order 
to build the minimal decomposition for every state has 
been given too. The classification of the pure three-qubit 
states in terms of their entanglement properties can be 
done following alternative criteria to the one described 
here, which is based on the features observed acting with 
the group of local unitary transformations. A possible 
approach is to classify the states looking for their prob- 
abilistic conversions under local operations and classical 
comunication (LOCC) for the single-copy case (see | [l^ ] 
and also ||lj,|l^) or in the asymptotic regime |1^]. It 
would be expected that these classifications are a coarse- 
graining of the one presented in this work. In fact this is 
the case for the equivalences classes under LOCC given 
in p. 

Finally it has been indicated how to extend decompo- 



sition (11) to systems of more parties. A simple counting 
of parameters shows that at least 2^ — N product states 
built from local bases are needed in order to specify a 
generic A^-qubit pure state, and for four qubits we suc- 
ceeded to find a procedure that makes zero four of the 
coordinates tijki- The decomposition (nSh allows for a 
simpler generalization. However, in all the cases some 
difficulties arise, related to the uniqueness of the decom- 
positions, and it is not clear whether these generalized 
Schmidt decompositions are quite useful for composite 
systems of more than three qubits. 



APPENDIX A: REAL STATES 

In this appendix we will show that, given a pure three- 
qubit state \ip) e C^(8)C^(8)C^, this state is real, i.e. there 
exists a product basis for which all coefficients are real, 
if and only if \fT1J2Jl — ^ or Aj = 0. 

Consider the case of a state | V') = X^i 7 k ^ijklijk) where 
all the tijk are real. Now we will follow the procedure de- 
scribed by the equations (R-foh that gives us the decom- 
position (11]). Since the initial coordinates are real, from 
(0) a second degree equation in x with real coefficients is 
obtained, and this implies that the two solutions, x and 
X, satisfy that either they are both real or x = x*. In 
the first case, the calculation of the decomposition can 
be performed using orthogonal matrices, and since the 



initial coordinates were real, we will obtain a real de- 
composition, i.e. ip = 0,Tr, which is equivalent to (p2). 
For the second case, since x = x* , ti{T^^T^) = tr(to'''t^), 
and then fiQ = i^o and ( p3|) is satisfied. 

Now, the inverse has to be proved. For the first case 
it is clear that all the states verifying ( p^ ) take real co- 
ordinates when they are expressed in the basis used in 
decomposition ([ll|). For the second case the proof is not 
so trivial. 

Consider a generic state, |(/)), having Aj equal to zero. 
The parametrization of this family of states is simplified 
using (p^), so let us first mention some facts about this 
decomposition. As it has been shown, any state with 
J4 7^ can be written as (09) pJlS] where 

1 






Ji 



/3 = -^ VAi2A*3+M4(A*4+M2+M3) 

M 
,5 = arg(AiA4e^'^-A2A3), 



(39) 



and, up to unitary transformations. 



|0> 



m 



coa-f, 
sin 7j 



i = l,2,3. (40) 



It can be proved that when Aj = the coefficients a and 
(3 are equal and then the states to be studied are 



a{\000) + e'^\ipiip2^3)) 



(41) 



Recall that for these states the complex conjugate is in 
the same orbit as the original one, and this means that 



ijk 



J2^^<^^%^kctabc, 



(42) 



where tijk are the coordinates in some product basis and 
and vf^ are the elements of the local unitary 



"ja; ""jb 



matrices, V^ in A, V^ in B and V^ in C, connecting 
the two states. From (^l]) it follows that these unitary 
operators are 



V = e- 



iS' / c 



(43) 



where Ci = cos7i, Si = sinji and 6' = -^ (actually, the 
phase factors in the matrices V^ can be given by arbi- 
trary angles Si satisfying the constraint ^^ 6i = S, but 
we choose these angles for simplicity). 

Now we would like to find a product basis for which 
all the coefficients are real, i.e. 



^ijk — Z^'^ia^jb'^kc^^abc — tijki 



Hjk. (44) 

and from this condition and using (l42), we have 

y* = {y/'^fwK (45) 



The explicit form of each V^ , ( [1 3D , as a product of a phase 
factor and a real and symmetric matrix allows to write 
them as 



a(|000)±|(^i(p2^3)). 



(51) 



V 



e-'^'{0' 



D'O' 



(46) 



where O* are orthogonal matrices and Z?* are diagonal 
matrices with entries ±1. The change of basis we are 
looking for then is given by 



W = (D^^C ^ e 



-iS" 



Ci 

-isi 



ici 



(47) 



where Ci = cos ^ , s,: 



sin ^ and 5" = 4 . The new 

2 6 



coordinates obtained applying these local change of basis 
are, up to normalization, 



ioOO = £lC253COS- 


*ooi — 


-CiC2S3sin- 


.... s 

tow = -ciS2C3Sin- 


^011 — 


6 
-C1S2S3COS- 


.... s 
iioo = -S1C2C3 sm - 


'^101 — 


S 
-S1C2S3COS- 


*110 = -S1S2C3COS- 


t'lii = 


S 
siS2S3Sin- 


is ends the proof. 







(48) 



APPENDIX B: TYPE 4D 

In this section we prove that a thrcc-qubit pure state 
\ip) can be written as 



1^) = ^ilOOl) + ^alOlO) + Z3|100) + ;4|111) 



(49) 



and 



with real and positive coefficients, if and only if (|3S 
( |33| ) are verified. 

Starting from ( |49|) we can apply the procedure given by 
(Kr|lO|) to obtain (|ll|). It can be seen that all the unitary 
matrices needed for the determination of this decompo- 
sition are real, i.e. they are orthogonal, and since the 
original coefficients {li} were also real, we will obtain a 
real decomposition with (p2h. Moreover, it can also be 
proved that the two matrices obtained after (@), Tg and 
Tq, corresponding to the two solutions of this equation, 
X and X, verify 



tr((r^)tT^)=tr((T^)tr^). 



(50) 



This condition implies that /io = p-o, and using (p3) we 
have also (p3|). 

Now we prove the inverse. Consider a state \(j)) satis- 
fying (B3) and (p3). Because of the latter condition, the 
state allows for a decomposition as (pl|). Moreover, since 
(B3) is also satisfied, we have (/s = 0,7r in (pi]), and this 
implies, using (p9|), that S — 0,7r. The generic expres- 
sion for a state satisfying both the conditions can be now 
given. 



If we perform the local change of bases described by ( |47| ) 
it can be seen, using (Ba) and the fact that S — 0,t:, that 
the state \(j>) is of type 4d. Indeed, the new coordinates 
are, after absorbing the phases and up to normalization. 



h = S1S2C3 I2 = S1C2S3 ^3 = C1S2S3 I4 



C1C2C3, 



(52) 



for 6 — 0, and 

h = C1C2S3 h = C1S2C3 ^3 = S1C2C3 I4 = S1S2S3, 

(53) 

for S = TT. Note that the local bases that appear in 
([49I ) are the ones that diagonalize the local density ma- 
trices. This gives the procedure to be applied in order to 
find the minimal decomposition without performing the 
maximization of (p4), which is generically a more difficult 
calculation. 
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FIG. 1. The figure depicts the three inequivalent sets of 
states given by (|), (§) and (|). 



Type 


Conditions 


States 


4a 


J4=0,VJlJ2J3 = f 


|000),|100),|101),|110) 


4b 


J2 = J5 = 


|000),|100),|110),|111) 


4c 


Ji J4 -f J1J2 + Ji J3 + J-iJ-i = 


|000),|101),|110),|111) 


4d 


Aj^0,VJlJ2J3^^ 


|001),|010),|100),|111) 


3a 


J1J2 + J1J3 + J2J3 = 
v'JlJ2J3 = ^,J4 = 


1000), 1101), 1110) 


3b 


Ji = J2 = Js = 


1000), 1110), 1111) 


2a 


All Ji — apart from Ji 


1000), 1011) 


2b 


All Ji = apart from J4 


1000), 1111) 


1 


J^^0 


1000) 



TABLE I. Classification of three-quantum-bit states. For 
the types of states denoted by 4b, 4c, 3b and 2a, there exist 
analogous condition interchanging the roles of the invariants 
Ji, J2, J3, and consequently the product states used in the 
minimal decomposition. 



